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The low-energy limits of models with disorder are frequently described by sigma models. In two 
dimensions, most sigma models admit either a Wess-Zumino-Witten term or a theta term. When 
such a term is present the model can have a stable critical point with gapless excitations. We 
describe how such a critical point appears, in particular in two-dimensional superconductors with 
disorder. The presence of such terms is required by the underlying (anomalous) symmetries of the 
original electron model. This indicates that the usual symmetry classes of disordered systems in 
two dimensions can be further refined. Conversely, our results also indicate that models previously 
thought to be in different universality classes are in fact the same once the appropriate extra terms 
are included. 



Understanding the effects of disorder is one of the cen- 
tral problems of condensed matter physics. It is natural 
to focus upon the phase transitions these systems un- 
dergo. In this paper, we will describe several ways phase 
transitions and the corresponding gapless modes arise in 
disordered models in two dimensions. 

The low-energy degrees of freedom for non-interacting 
fermions with disorder can frequently be described by 
a sigma model in the replica formulation (![]. A G/H 
sigma model arises when the vacuum manifold of a field 
is H, a subgroup of the global symmetry group G. The 
low-energy modes of this model then take values in the 
coset G/H. Such sigma models in two dimensions do not 
generically have gapless degrees of freedom. For finite 
number of replicas, this is a result of the Mermin-Wagner- 
Coleman theorem. The sigma models are usually gapped 
because the space G/H is curved. Expanding the action 
in powers of the curvature gives a mass scale: this is the 
well-known story of asymptotic freedom and dimensional 
transmutation. At specific values of N, the curvature and 
hence the beta function can vanish Q , but in this paper 
we will be concerned with the general situation. 

Despite the fact that sigma models are usually gapped, 
there are many non-trivial fixed points in two dimensions. 
One can add to the sigma model action terms which dra- 
matically change the low-energy physics. One famous 
example is the topological theta term. For example, the 
model with G = 0(3) and H = 0(2) has a gap || unless 
one adds a theta term to the action. This 0(3)-invariant 
term is a total derivative, so it depends only on the fields 
at spacetime infinity and has no effect on perturbation 
theory around the unstable Gaussian fixed point. Ex- 
plicitly, the theta term takes the form inO, where n is 
an integer counting how many times a field configura- 
tion winds around the 0(3)/0(2) two-sphere. As has 
been discussed in the context of the Heisenberg spin- 
chain, at 9 — 7T there is a non-trivial low-energy fixed 
point [[H|. It is unstable when 9 is not tt, so tuning 9 
through tt results in a phase transition. This behavior 



inspired the well-known proposal that the Hall plateau 
transition is described by N — > replica limit of the 
U(2N)/U{N) x U(N) sigma model with a 9 term §. 

The theta angle need not be a tunable parameter. If 
there is a symmetry of the model (like time-reversal or 
CP invariance) under which the theta term is odd, then 
9 is fixed at zero or tt. Moreover, in certain sigma mod- 
els such as SU (N) / SO(N), the winding number can only 
be or 1, so 9 can be only zero or tt. As we will discuss 
below, in such situations 9 is determined by the underly- 
ing disordered model. When 9 = tt, the SU{N)/SO(N) 
model has a stable low-energy fixed points: it flows to 
SU(N)i Q. Thus the fixed value of 9 determines whether 
the model is gapless or not. 

There is another way to modify some two-dimensional 
sigma models to obtain a stable low-energy fixed point. 
This is to add a Wess-Zumino-Witten (WZW) term to 
the action 0. In four-dimensional particle physics, the 
WZW term for example describes the tt° — > 2j decay 
in low-energy effective field theory JtJ. Such terms arise 
when the fields take values in a group. Consider a theory 
with global symmetry H x H, where H is any simple Lie 
group (e.g. SU(N),SO(N), or Sp(2N)). If the vacuum 
manifold of the theory is chosen such that the modes in 
the diagonal "vector" subgroup Hy become gapped, the 
low-energy fields, h, take values in the remaining space, 
H x H/Hy- This space is isomorphic to H itself, and 
is called the axial part of the original symmetry group. 
For example, for SU(N), h would be an N x N unitary 
matrix with determinant one. 

Without the WZW term, a sigma model where the 
fields take values in a simple Lie group is called a princi- 
pal chiral model. In terms of the coupling, g (the inverse 
stiffness), its action is 

Spcm = ~ J dxdy Tt [d^h d^hT 1 ] . (1) 

We assume that the fields all fall to a constant at spa- 
tial infinity, so that the two-dimensional spacetime can 
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be treated as a sphere. To write the WZW term, one 
needs to extend the fields h(x, y) on the sphere to fields 
h(x,y,z) on a ball which has the original sphere as a 
boundary. The fields inside the ball are defined so that 
h at the center is the identity matrix, while h on the 
boundary is the original h(x,y). It is possible to find a 
continuous deformation of h(x, y) to the identity because 
TaCff) = f° r an Y simple Lie group. Then the WZW 
term is kT, where 

r= ^T I dxd V dzTT [{h- 1 d i h){h- 1 d :j h)(hr' i -d k h)\ . 

ZQ7T J 

(2) 

The coefficient k is known as the level, and for compact 
groups must be an integer because the different possible 
extensions of h(x, y) to the ball yield a possible ambiguity 
of 2-7T times an integer in T. Unlike the 9 term, the WZW 
term does change the equations of motion, but only by 
terms involving h(x,y): the variation of the integrand is 
a total derivative in z. 

The sigma model with WZW term has a stable fixed 
point at g = 16n/k ||], so the model is critical and the 
quasiparticles are gapless. The corresponding conformal 
field theory is known as the H k WZW model ||. The 
WZW term is invariant under discrete parity transforma- 
tions (e.g. x —> —x, y —> y) provided h is simultaneously 
transformed via h — > h~ l . In other words, the WZW 
term in a parity-invariant theory involves pseudoscalars. 

The WZW and theta terms are closely related. For ex- 
ample, if one adds a term \{Trh) 2 to the SU(2) k WZW 
action, at A large the low-energy modes live on the sphere, 
i.e. the SU(2)/0(2) sigma model. One can check explic- 
itly that the WZW term then turns in to the theta term 
with 9 — kir [Q. Moreover, in sigma models with only 
winding number or 1 , the topological term can be writ- 
ten in the form of the WZW term @ |j. 

Given the above considerations, it is important to un- 
derstand how a WZW term arises in the replica formula- 
tion of two-dimensional models with disorder. We illus- 
trate the appearance first by considering a simple theory 
where such a term appears, related to models studied in 
P,[lO[ . We start with a system of spin-polarized fermions 
with a triplet p-wave type pairing: 



Hq = '^2 e k c f k c k + (A k c' k &_ k + h.c), 

k 



(3) 



where A k ~ ^k x . We study disorder weak enough to 
maintain some notion of the Fermi surface. The low en- 
ergy excitations of the fermions are then found about two 
nodes positioned on the fc y -axis at K± = (0, ±A). We 
then linearize the theory about these nodes via tK ± +q — 
q y vp and c ~ c\ exp(iK + x) + C2 exp(iK-x). Introducing 
the spinor ip^ = (cj, 0%), we can write H as 



H 



d 2 x il)\iv F T z d y + iv A T x d x )ijj. 



(4) 



The Pauli matrices t, act in the particle-hole space of the 
spinors. H is precisely the Hamiltonian of a single Dirac 
fermion in 2 + 1 dimensions. To obtain the sigma model 
description of this Hamiltonian, we fix the energy, u>, at 
which we work, and describe the theory in terms of an 
action of a two-dimensional Euclidean field theory: 



So 



d 2 x ip'(ivFT x d y + ivAT x d x - iu>Tz)ip. (5) 



With this action we can compute correlators of the form 
(l/(H-u)). 
We include on-site disorder by adding a term, 



H, 



disorder 



t(c\ Cl 



c\c 2 ) 



(6) 



to the Hamiltonian. Here t is a random variable with 
variance (t(x)t(y)) = (2u)~ 1 5(x — y). To compute cor- 
relators with disorder we replicate the theory, if> — > ip k . 
Disorder is then easily averaged over leaving a set of quar- 
tic terms, ^disorder = ~l i ('V k T z ^k)(i>\T z '4)i). Reorganiz- 



ing the fields via = (ip R , W L ) — ^t z exp(i7rr :c /4) and 



tp = (ipRjipi,) = exp(— inT x / , gives the action (at 
uj = 0) to be: 



S 



1 < 



d'xi^l (v F d y -v A r z d x )i, k -- (V£ Vfc ) ) • (7) 



This theory is invariant under the group U(N) l X U(N)r 
transforming ip -* |((1 + t z )U l + (1 - T z )U R )tp. Adding 
disorder to the Cooper-pairing term results in another 
four-fermion term, but preserves this symmetry. As long 
as the symmetry structure is unchanged, the low-energy 
sigma model should be the same. 

When to = 0, this action is equivalent to a massless 
Dirac fermion in a random magnetic field; the chiral 
symmetry prevents a mass term from appearing. This 
model can be solved by bosonization p2| . Free fermions 
are equivalent to the WZW model SU(N) 1 x U(l) 
The four-fermion term is then simply expressed as 
J(Jl + Jr) 2 , where Jl and Jr are the U(l) currents, 

Jl = {^l^l) and J R = Wr^r)- These ^(1) currents 
can be expressed in terms of a free boson with ji, = id^ 
and jn = —idu<j). The four-fermion coupling then merely 
determines the boson radius. Therefore, the model re- 
duces to pure SU(N)i together with a decoupled free 
boson, and hence is critical. 

To see directly how the WZW term appears, it is use- 
ful to rederive this result from the sigma model approach 
pH . We assume there is an energy scale where some 
fermion bilinear gets an expectation value. Formally 
speaking, we introduce a Hubbard-Stratonovich matrix 
field, M k i, to factor the four-fermion term: 



S = So - j d 2 x(-Tr(M 2 ) + i^lMu^)), (8) 

where M is hermitian. Under U(N) L x U(N) R , M -> 
UMU^ . The WZW term appears in low-energy effective 
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theories describing saddle points where M is off-diagonal, 
e.g. M LL = M RR = 0, but M LR = M RL oc I, the iden- 
tity. The diagonal subgroup U{N)v leaves this saddle 
point invariant, and so the low-energy modes T = M R l 
take values in U{N) L x U(N) R /U(N) V « U(N). Note 
that under parity, «-> T. Focusing solely upon these 
modes allows us to write 

3 = So - J rf 2 x(^Tr(TTt) + i[^ R T^ L + ^ L T^ R )). 

(9) 

Integrating out the fermions leaves an effective action for 
the bosonic field T, which can be expanded in powers of 
the momentum over the expectation value of T. At low 
energy, one can safely neglect four-derivative terms and 
higher, and one easily obtains the kinetic term (Q) for 
T. However, a less-obvious term also appears. The case 
of interest was treated in Q]. To perform a consistent 
low-energy expansion, one must change field variables. 
This results in a Jacobian in the path integral [jl5| , which 
is precisely the WZW term. Writing T = he^, where 
det(h) = 1, we obtain here the WZW term (||) of the 
SU(N) field h, with level k = I. The extra U(l) field 4> 
in the sigma model is a theory of a decoupled boson. 

Thus in the sigma model approach we recover the crit- 
ical line described by SU (N)i xU(l) conformal field the- 
ory in the limit N — > 0. This critical line is a part of the 
critical space of the Gade hopping model; the full space 
is larger because the SU (N) coupling does not flow when 
N = [p|ptl- Moreover, this same CFT also applies to 
a time-reversal symmetric version of the hopping model, 
which is described by a U (N)fO(N) sigma model. When 
9 = 7r, this sigma model flows to the same SU (N)i x U (1) 
conformal field theory ||. We can see how the two 
are related by perturbing the U(N) WZW model ac- 
tion by a term which gives a mass to the modes out- 
side a U(N)/0(N) subspace (the subspace corresponds 
to symmetric U(N) matrices). As with the N = 2 case 
mentioned above, the low energy limit of the perturbed 
U(N) WZW model reduces to the U(N)/0(N) sigma 
model with 9 = kit. Thus the Gade model with and 
without time reversal invariance provides a concrete re- 
alization of the equivalence of theta and WZW terms in 
disordered systems. 

In no sense does the WZW term arise in the above 
models as a result of fine tuning: it must appear when 
there is a chiral anomaly. In the model treated above, 
the fermions have a U(N)l x U(N) r symmetry. As is 
well known, chiral symmetries involving fermions are fre- 
quently anomalous, so that the Noether currents do not 
all remain conserved in the quantum theory. For mass- 
less fermions in 1 + 1 dimensions, this was shown in detail 
in |l7|fi|. The WZW term is the effect of the anomaly on 
the low-energy theory. Even though the fermions effec- 
tively become massive when T gets an expectation value, 



their presence still has an effect on the low-energy the- 
ory, even if this mass is arbitrarily large. This violation 
of decoupling happens because the chiral anomaly must 
be present in the low-energy theory. In other words, the 
anomaly coefficient does not renormalize. This follows 
from an argument known as 't Hooft anomaly matching 
]l6[ . One imagines weakly gauging the anomalous sym- 
metry (in our case the axial SU (N)). It is not possible to 
gauge an anomalous symmetry in a renormalizable the- 
ory, but one can add otherwise non-interacting massless 
chiral fermions to cancel the anomaly. Adding these spec- 
tator fermions ensures that the appropriate Ward iden- 
tities are obeyed and the symmetry can be gauged. In 
the low-energy effective theory, the Ward identities must 
still be obeyed and the theory must remain anomaly-free. 
Because the massless spectator fermions are still present 
in the low-energy theory, there must be a term in the 
low-energy action which cancels the anomaly from the 
spectators. This is the WZW term. 

We now turn to a case of considerable interest, that 
of d x 2_ y 2 superconductors in the presence of disorder 
Jl8| , p~9| | . We now study fermions with spin symmetry, 
so that c, c' carry a node index, a spin index, and a 
replica index. We label the nodes (1, ±) for the pair 
near the wavevectors fc = ±(ir/2, ir/2) and (2, ±) for the 
pair near k = ±(— 7r/2, 7t/2). It is convenient to group 
particles and holes together. Thus define four doublets, 
i4ux = (c 3 aa ,i(cr v s x c : >' l ) aa ), a = ±,a =t / |, for each 
pair of nodes j — 1,2. The Pauli matrices a a and s a act 
upon the spin/± node indices respectively. If we rotate 
the system in the ab plane by 45°, the linearized Hamil- 
tonian for a superconductor with time-reversal and spin 
symmetry is 

H = hp\v F s z T z d x + v A s z T x dy)tp 1 + (x, 1 y, 2), (10) 

where the Pauli matrix r acts in particle/hole space 

and ip = (s x <r y t v ipi ) T . After adding disorder, the 
replicated Hamiltonian is invariant under the group 
Sp{2N) L x Sp(2N) R @. It sends V' -> Uip\ where 
U = i[f/i(l+r") + U r {1-tV)], with U L and U R 
each elements of Sp(2N) acting only upon spin and 
replica indices. An element Q of Sp(2N) is an in- 
vertible 2N x 2N real matrix obeying Q T a v Q = a y . 
With Gaussian on-site disorder and uj = 0, the low- 
energy modes of the corresponding action take values in 
Sp(2N) L x Sp{2N) R /Sp{2N) v . The low-energy action 
thus includes the term @) with h in Sp(2N) Q. 

This model also has a WZW term. We first consider 
the case where disorder does not couple the pairs of 
nodes. Once disorder has been averaged over and the 
quartic terms have been factorized, the fermions interact 
with hermitian matrix fields, and M( 2 \ via terms 

of the form 

^iM^l+^ a M^j, (11) 
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as before. At a saddle point with non-zero values of the 
off-diagonal elements, the argument of jL4j shows that we 
will again obtain level 1 WZW theories for M W and for 
M&\ in this case Sp(2N) 1 . 

We make this explicit by showing how the spin degrees 
of freedom expand the previous U (A) symmetry to that 
of Sp(2N). Sp(2N) has a U(N) subgroup consisting of 
orthogonal real matrices of the form 

A B 
-B A 

where we put the up spins in the first N components of 
ip, and the down spins in the second half. Thus the U(N) 
subgroup acts on the combinations, ip^ ± itpi, with the 
matrices A =p iB. In group-theoretic language, the 2N 
dimensional representation of Sp(2N) decomposes into 
the A + A representation of U (A) . We can examine the 
WZW model describing the low-energy excitations which 
lie in the U (A) subspace. It must be exactly the same as 
that treated above, with the exception that there are now 
two species of "spinless" fermions ± itpi . This results 
in two WZW terms, one with h and one with h*. As the 
WZW term is invariant under complex conjugation, the 
two contributions are identical, and so add. This the- 
ory is a U(N)2 WZW theory (as was also derived using 
bosonization Jl8| ; the corresponding result in the super- 
symmetric formulation was derived in p4j ) . The level of 
the WZW term in the full Sp(2N) theory is determined 
by the embedding: when an Hi theory is embedded in 
Gfe, the levels obey kr = /, where r is a group-theory fac- 
tor called the index of the embedding of H into G. For 
the embedding of U(N) into Sp(2N), r = 2. Thus we 
have shown that the low-energy theory for ([ll]) is given 
by an Sp(2N)i WZW model for each pair of nodes. 

When a WZW term is present, the quasi-particles of 
the spin/time-reversal invariant superconductor need not 
be localized. Specifically, the coupling constant of the 
Sp(2N)k theory, inversely related to the spin conduc- 
tance, now flows to a fixed value (provided its bare value 
is sufficiently small), as opposed to becoming arbitrar- 
ily large. Thus with a WZW term such models are no 
longer spin insulators. The Sp(2N)k WZW models seem 
to be well behaved as A — > 0. The dimensions of the 
fundamental operators are given by m(2 — m)/(2k + 2), 
where m is a positive integer. To compute the density of 
states at a finite energy, u>, we add a term of the form 
uTrQi + W) to the action. The density of the states is 
then given by p(uj) ~ (Tr(h + tf)). This field tr h has 
dimension A = l/(2fe+2), and a simple scaling argument 
yields p{u) ~ u A /( 2 - A ) = w 1 /^ 3 ), This result agrees 
with a previous computation [ fl8[ when k = 1. 

When the two pairs of nodes are not coupled by the 
disorder, the low-energy theory is therefore described by 
two Sp(2N)i WZW models, which have symmetry cur- 
rents, Jli,2, Jri,2- Once the pairs of nodes are coupled, 
there are three possibilities. The first is that the two 



Sp(2N)l x Sp(2N)n symmetries are preserved, so the 
low-energy theory remains two Sp(2N)i WZW models. 
The second is that only simultaneous chiral transforma- 
tions on the two pairs remain a symmetry. In the decou- 
pled theory, this symmetry is generated by J^i + Jn 
and Jrj + Jr 2 . It has level 2 (in the language used 
above, the index of the embedding is 2). Once the nodes 
are coupled, the generators may be deformed away from 
Jli + Jl2 and Jri + Jm, but the level cannot change. 
This occurs for the same reasons that anomalies remain 
in the low-energy theory - the anomaly coefficient cannot 
renormalize. Thus the low-energy behavior of this model 
is described by an Sp{2N)2 WZW model. If one modi- 
fies the Hamiltonian so that the two pairs of nodes are 
no longer identical, then the result is the same as long as 
the chiral symmetry is not explicitly broken: the chiral 
anomaly does not change under perturbation. 

The third possibility is that the anomalies cancel be- 
tween the two theories. This is in fact what happens 
in d x 2_ y 2 superconductors (2l). Consider two coupled 
WZW models with fields hi and hi and action 

SpcMih^+Tih^ + SpcM^+T^ + XTrih^ + h.c). 

The global symmetries of the decoupled model are hi — > 
ULihiUfti. In the coupled model this global symmetry is 
broken to h\ — > Ui,ihiU' L2 , hi — > UL2hil/j jl , appropriate 
to the underlying lattice symmetry of the d-wave super- 
conductor. These symmetries are anomaly free; the axial 
anomalies of the two nodes cancel upon coupling. This 
can be seen roughly in that the coupling A is relevant; in 
the strong coupling limit h\ — h\ (i.e. r(/ii) = — T{h<i)) 
and so the WZW terms cancel. This omits terms that 
might arise in integrating out massive degrees of freedom. 
This possibility cannot be ignored, because we know this 
is precisely how a WZW term arises in . However, this 
argument can be made more precise in the case when h\ 
and hi are elements of O(N) instead of Sp(2N). Indeed, 
then one can fermionize one of the two WZW models, 
writing (/i 2 ) a h = Tpbipa, with the resulting action 

SpcmQix) + T(hi) + + \@hnl> + h.c). 

For A large, we can integrate out these fermions accord- 
ing to |14j as before. In doing so, a WZW term for hi is 
induced which precisely cancels the original WZW term 
for h\. If we had instead coupled the theories via a term 
Tv(hih\ + h.c), the two terms would have added, and 
we would have recovered the level two theory discussed 
above. This cancellation occurs even if the two theories 
have differing Fermi velocities. Including Fermi velocities 
directly in the action (easily done through dimensional 
analysis) shows that the terms Spcm are affected while 
the WZW terms T are not. Thus the cancellation be- 
tween the two WZW terms can proceed as above. 

A d-wave superconductor with spin rotational but bro- 
ken time reversal invariance is realized when the gap 
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wave function is A 



id r 



This breaks the 



Sp(2N) x Sp(2N) symmetry explicitly, and the appropri- 
ate sigma model in the replica formulation for this theory 
is Sp(2N)/U(N) with a theta term Because the 

time-reversal invariance is broken, the coefficient of 9 is 
not constrained by symmetry, but at 9 = it one expects a 
critical point in the same fashion as H. It is conjectured 
in Q that at 9 — n, this model flows to the Sp(2N)i con- 
formal field theory. By utilizing the supergroup formula- 
tion of the corresponding network model (spin quantum 
Hall effect), it was shown that the critical point point of 
the disordered model is equivalent to classical percolation 
p|. Thus the replica limit N -> of Sp(2N) 1 should be 
equivalent to the critical point in the spin quantum Hall 
effect. The exponents computed above do indeed agree 
with those computed in |23|| , for example p{E) cx E 1 / 7 . 
We again see a feature that appeared in the U(N) class 
of problems: related models seemingly in two different 
universality classes are in the same one, once the WZW 
terms and 9 terms are specified. 

We have shown how stable critical points can arise in a 
number of replica sigma models. This shows that a num- 
ber of disordered symmetry classes in two dimensions [ |24| ] 
allow further refinement. In particular, in models with a 
chiral symmetry, one must compute the anomaly in the 
original underlying short-distance replica (or supergroup) 
model to determine the WZW term in the low-energy 
model. Likewise, in sigma models with Z2 winding num- 
bers, one must determine whether 9 — or tt. 
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